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Shahar Hod 

The Ruppin Academic Center, Emeq Hefer 40250, Israel 

and 

The Hadassah Institute, Jerusalem 91010, Israel 
(Dated: November 24, 2008) 

We study analytically the relaxation phase of perturbed, rapidly rotating black holes. In 
particular, we derive a simple formula for the fundamental quasinormal resonances of near-extremal 
Kerr black holes. The formula is expressed in terms of the black-hole physical parameters: uo = 
m£l — i2irTBH(n+ |), where Tbh and fl are the temperature and angular velocity of the black hole, 
and m is the azimuthal harmonic index of a co-rotating equatorial mode. This formula implies that 
the relaxation period r ~ 1/Qu) of the black hole becomes extremely long as the extremal limit 
Tbh — » is approached. The analytically derived formula is shown to agree with direct numerical 
computations of the black-hole resonances. We use our results to demonstrate analytically the fact 
that near-extremal Kerr black holes saturate the recently proposed universal relaxation bound. 
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The radiative perturbations of a complete gravitational collapse decay with time leaving behind a 'bald' black hole. 
This is the essence of the no-hair conjecture introduced by Wheeler [l| more than thirty years ago. It asserts that 



perturbation fields left outside the collapsing star would cither be radiated away to infinity, or be swallowed by the 
newly born black hole. 

According to the uniqueness theorems @, 0, 0, H, @] , the metric outside the black hole should relax into a Kerr- 
i— i' Newman spacetime, characterized solely by the black-hole mass, charge, and angular momentum. This relaxation 
phase in the dynamics of perturbed black holes is characterized by 'quasinormal ringing', damped oscillations with a 
discrete spectrum (see e.g. [7| for a detailed review). At late times, all perturbations are radiated away in a manner 
reminiscent of the last pure dying tones of a ringing bell [1, 0, [l(J EH • 

The black hole quasinormal modes (QNMs) correspond to solutions of the perturbations equations with the physical 
boundary conditions of purely outgoing waves at spatial infinity and purely ingoing waves crossing the event horizon 
Such boundary conditions single out a discrete set of black-hole resonances {w„} (assuming a time dependence of 
the form e~ luJt ). In analogy with standard scattering theory, the QNMs can be regarded as the scattering resonances 
of the black-hole spacetime. They thus correspond to poles of the transmission and reflection amplitudes of a standard 
QO 1 scattering problem in a black-hole spacetime. 

The characteristic quasinormal frequencies are complex. This reflects the fact that the perturbations fields decay 
with time, in accord with the spirit of the no-hair conjecture. It turns out that there exist an infinite number of 
quasinormal modes, characterizing oscillations with decreasing relaxation times (increasing imaginary part) [131 ]. The 
mode with the smallest imaginary part (known as the fundamental mode) determines the characteristic dynamical 
timescale r for generic perturbations to decay. 

Quasinormal resonances are expected to play a prominent role in gravitational radiation emitted by a variety of 
astrophysical scenarios involving black holes. Being the characteristic 'sound' of the black hole itself, these free 
oscillations are of great importance from the astrophysical point of view. They allow a direct way of identifying the 
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spacetime parameters, especially the mass and angular momentum of the black hole. This has motivated a flurry of 
research during the last four decades aiming to compute the quasinormal mode spectrum of various types of black- hole 
spacetimes Q- 

It is well known that realistic stellar objects generally rotate about their axis, and are therefore not spherical. Thus, 
an astrophysically realistic model of wave dynamics in black-hole spacetimes must involve a non-spherical background 
geometry with angular momentum. In this work we determine analytically the fundamental (least-damped) resonant 
frequencies of such rap idly-rotating Kerr black holes. (For a recent progress in the study of the highly-damped 
resonances, see Il5|.) The spectrum of quasinormal resonances can be studied analytically in the near-extremal 
limit (M 2 — a 2 ) 1 / 2 -C a < M, where M and a are the mass and angular momentum per unit mass of the black hole, 
respectively. 

Before going on we would like to summarize what is already known about these fundamental, slowly-damped Kerr 
QNMs: 

• In the extremal limit (a — > M) one finds ~ m/2M, where m > is the azimuthal harmonic index of a 
perturbation field co-rotating with the black hole. This fact has been found in many numerical computations 
(see e.g., 13, 16, 17, HH) and is well understood analytically H, 2ll 21 \- 
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• Numerical computations ID, 3, 12, El have indicated that, for co-rotating modes — > in the extremal 
limit. This implies that these modes are long lived. Detweiler presented a semi-analytical formula for 
the long-lived black-hole resonances in the extremal limit. It is important to emphasize that Detweiler's result 
assumes that ^slu/Tbh >> 1 and is not appropriate for the most long-lived modes in the non-extremal case, 

where it might be that 9w ~ Tbh , where Tbh = ^MlM+iM^-a 2 ) 1 / 2 ] ^ s * ne Bekenstein-Hawking temperature of 
the black hole. 

• Motivated by the recently proposed time-temperature universal relaxation bound (TTT) [22I ] , we have recently 
re-examined the role played by the imaginary parts of the black-hole quasinormal resonances. It has been 
observed [22I I2H ] that the numerically computed equatorial resonances of near-extremal Kerr black holes are 
well-approximated by the simple analytical relation lo — mfl — Vl^Tsnin + 1/2), where n = 0, 1,2, .... The 
successfulness of this conjectured formula is demonstrated in Table Q] The numerical results therefore imply 

O (Tbh) in the near-extremal limit. A similar relation which is valid only for the fundamental mode was 



obtained in 24| for equatorial modes. Table HT1 demonstrates the fact that the predictions of the formula improve 



as the extremal limit is approached. 



n 




Qlo (numerical) 


Quj (analytical) 





0.99324 


0.00341 


0.00348 


1 


0.99322 


0.01020 


0.01045 


2 


0.99321 


0.01699 


0.01743 


3 


0.99320 


0.02385 


0.02440 


4 


0.99317 


0.03067 


0.03137 


5 


0.99313 


0.03749 


0.03834 



TABLE I: Quasinormal resonances of a near-extremal Kerr black hole with a/M = 0.9999. The data shown is for the equatorial 
mode I = m — 2, see also [l7]]- The proposed analytical formula is lu — mfl — ¥1-kTbh(ti + 1/2). (Here mQ ~ 0.9859). The 
agreement between the numerical data and the proposed formula is of ~ 2%. 



a/M 




^SLOana / ^'(-^num 


0.9 


0.933 


1.170 


0.96 


0.977 


1.106 


0.9999 


0.993 


1.022 



TABLE II: The ratio between the relation lu = mil — %-kTbh for the fundamental equatorial black-hole resonance, and the 
numerically computed value. The data shown is for the mode I = m = 2. The agreement between the numerical data and the 
analytical formula improves as the black hole approaches its extremal limit. 



In order to determine the black-hole resonances we shall analyze the scattering of massless waves in the Kerr 
spacetime. The dynamics of a perturbation field ^ in the rotating Kerr spacetime is governed by the Teukolsky 
equation [251 ]. One may decompose the field as (we use natural units in which G = c = h = 1) 



S s i m {6\au>)il) s i m (r)e~ 



(1) 



where (t,r,9,<f>) are the Boyer-Lindquist coordinates, lu is the (conserved) frequency of the mode, / is the spheroidal 
harmonic index, and m is the azimuthal harmonic index with —l<m<l. The parameter s is called the spin weight 
of the field, and is given by s — ±2 for gravitational perturbations, s = ±1 for electromagnetic perturbations, s = ±i 
for massless neutrino perturbations, and s = for scalar perturbations. (We shall henceforth omit the indices s,l,m 
for brevity.) With the decomposition |T]), %j) and S obey radial and angular equations, both of confluent Heun type 
[2^, [27]], coupled by a separation constant A(au). 

The angular functions S(9; au) are the spin-weighted spheroidal harmonics which are solutions of the angular 



equation 25, 27]] 
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The angular functions are required to be regular at the poles 9 = and 9 = tt. These boundary conditions pick out a 
discrete set of eigenvalues Ai labeled by an integer I. [In the atu <C 1 limit these angular functions become the familiar 
spin-weighted spherical harmonics with the corresponding angular eigenvalues A = 1(1 + 1) — s(s + 1) + O(a 2 oj 2 ).] 
The radial Teukolsky equation is given by 

. s cZ/ A , +1 AA\ \K 2 - 2is(r - M)K 2 , „ 4 „ . i, „ 

A s — ( A s+1 -M + j a 2 u 2 + 2mauo - A + Aisujr ij) = , 3 

or V dr J L A J 

where A = r 2 — 2Afr + a 2 and if = (r 2 + a 2 )w — am. The zeroes of A, r± = M ± (M 2 — a 2 ) 1 / 2 , are the black hole 
(event and inner) horizons. 

For the scattering problem one should impose physical boundary conditions of purely ingoing waves at the black- 
hole horizon and a mixture of both ingoing and outgoing waves at infinity (these correspond to incident and scattered 
waves, respectively). That is, 

( e -i*>v + K{uj)e lu y as r -► oo (y -> oo) ; 
\T{u)e~^ mn ^ &sr-^r+ (y -> -oo) , 

where the "tortoise" radial coordinate y is defined by dy = \(r 2 + a 2 )/A]dr. Here ft = 2A a [r+ is the angular velocity 
of the black-hole horizon. The coefficients T(uj) and 1Z(uj) are the transmission and reflection amplitudes for a wave 
incident from infinity. The discrete quasinormal frequencies are the scattering resonances of the black-hole spacetime. 
They thus correspond to poles of the transmission and reflection amplitudes. (The pole structure reflects the fact 
that the QNMs correspond to purely outgoing waves at spatial infinity.) These resonances determine the ringdown 
response of a black hole to outside perturbations. 

The transmission and reflection amplitudes satisfy the usual probability conservation equation |T(w)| 2 + |7?.(ct;)| 2 = 1. 
Teukolsky and Press [28| and also Starobinsky and Churilov [2^] have analyzed the black-hole scattering problem in 
the double limit a — ► M and ui — ► mfi. Detweiler [l9| then used that solution to determine the long-lived black-hole 
resonances in the extremal limit. Define 

a = ; t = M(lo — mVl) ; Co = ujr + . (5) 



Then the resonance condition obtained in [19( for a « 1 and r << 1 is: 

T(2id)T(l + 2i5)T(l/2 + s-2icj-iS)r(l/2-s-2iu>-i5) , , 0iK T(l/2 + 2iu + iS - Mr I a) 

' -2iua)™ - ' ■ — - (6) 



T(-2i5)T(l-2iS)T(l/2 + s-2iLu + iS)T{l/2-s-2iLu + i5) y ' r(l/2 + 2iCo - id - Air / 'a) 

where S 2 = Auj 2 — 1/4 — A — a 2 uj 2 + 2maoj. In the extremal limit Detweiler solved this equation with the assumption 
that t/<t — > oo. For near extremal black holes, it should be emphasized that the numerical data presented above is 
consistent with the limit t/cj — > const as a — > 0. Thus the quasinormal frequencies obtained in [19( do not include the 
most long lived, fundamental resonances of near-extremal black holes. We now derive analytically the fundamental 
resonances of rapidly rotating, near-extremal Kerr black holes. 

The left-hand-side of Eq. ([6]) has a well defined limit as a — > M and lj — > mfL We denote that limit by C. In 
the limit u) — » mfi, where u> is almost purely real, one finds from Eq. (|2|) that the separation constants {A} are also 
almost purely real. This in turn implies that the <5 2 's are almost purely real. For some modes, including most of the 
equatorial I = m modes (and also other modes with m close enough to /), 6 2 is found to be positive, which implies 
that in these cases 5 is almost purely real and positive [28|, |30(. For the rest of the modes one finds 5 2 < 0, which 
implies that 5 is almost purely imaginary with positive imaginary part (28l . |30| . 

If S is almost purely real and positive then one has (-i)- 2iS = e - 2i<51n (-«) = e -2i5in e -"/ 2 _ e -2iS^/2) = e -irS < L 
Here we have used the fact that S > 2 for all gravitational equatorial modes 28|, 3(3]. (In fact, this is also true for 



many other modes for which I~m>l.) If S is almost purely imaginary with a positive imaginary part then one 
has a~ 2lS — > in the near-extremal limit a — > 0. In both cases one therefore finds e = (— 2iLJcr)~ 2 ' 15 <C 1. Thus, a 



consistent solution of the resonance condition, Eq. ([6]), may be obtained if l/r(l/2 + 2icb ~ iS — Air/cr) = 0(e) [31 1. 
Suppose 

1/2 + 2iu> -iS- Air/a = -n + rje + 0(e 2 ) , (7) 
where n > is a non-negative integer and n is an unknown constant to be determined below. Then one has 

T(l/2 + 2iCo -iS- Air/a) ~ T(-n + T]e) ~ (-n) _1 r(-n + 1 + ne) ~ • •• ~ [(-l)"n!]- 1 r(?7e) , (8) 
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where we have used the relation T(z + 1) = zT(z) Next, using the series expansion l/T(z) = J2kLi c k zk with 

Ci = 1 [see Eq. (6.1.34) of [HI]], one obtains 

l/r(l/2 + 2iu) -id- Air/a) = (-l) n n\r]e + 0(e 2 ) . (9) 

Substituting this into Eq. © one finds r] = C/[(-l) n nW(-n + 2id)}. 

Finally, recalling that At /a = (ui — mQ) /2ttTbh we obtain from Eq. (JT]) the resonance condition 

(uj - mfl)/27rT B H = i[-n + r?e - 1/2 + i(S - to)] , (10) 

where we have substituted 2iu> ~ im for u> ~ mQ, ~ m/2M. The black- hole quasinormal resonances of equatorial 
Z = to > modes (and, in general, co-rotating modes for which 8 2 > 0) are therefore given by the leading-order 
formula 



uj = mQ, - i2irT BH (n + 1/2) + 0{T BH ) , (11) 

where n = 0, 1, 2, ... . One thus finds that $iu> — > mQ and 3cj — > in the near-extremal. These analytical findings are 
in accord with direct numerical computations, see Tables I and II [33| . 

The black-hole quasinormal resonances of non-equatorial I =/= to > modes (and, in general, co-rotating modes for 
which S 2 < 0) are given by the leading-order formula [see Eq. (TTU)) ] 

uj = mn- i2irT BH {n + 1/2 - iS) + 0(T BH ) , (12) 

where n = 0, 1, 2, .... It is worth pointing out that, in this case Qui > 2irTBH(n+l/2) (recall that 3(5 > 0). This implies 
that for these modes Qui approaches zero slower as compared to the equatorial ones for which Qui — 2irTBii(n+ 1/2). 
We have therefore established the fact that non-equatorial modes decay faster than the equatorial ones. 

In summary, we have studied analytically the quasinormal spectrum of rapidly-rotating Kerr black holes. The main 
results and their physical implications are as follows: 

(1) We have shown that the fundamental resonances can be expressed in terms of the black-hole physical parameters: 
the temperature Tbh, and the angular velocity fl of the horizon. 

(2) It was found that for all co-rotating modes (modes having to > 0) SRw — > mVL in the near-extremal limit. This 
conclusion is in agrement with, and generalizes, the I = to result obtained in (lij . 

(3) We find that, in the near-extremal limit Qui approaches zero linearly with the black-hole temperature. Namely, 
Qui = 0(Tbh)- This conclusion holds true for all modes co-rotating with the black hole. Thus, all co-rotating modes 
become long-lived as the black hole spins up. Moreover, it is realized that equatorial / = to modes (and in general, 
modes for which the quantity S is real) decay slower than other non-equatorial perturbations. 

(4) It is worth mentioning that a fundamental bound on the relaxation time r of a perturbed thermodynamical 
system has recently been suggested 22], t > h/itT, where T is the system's temperature. Taking cognizance of this 
relaxation bound, one deduces an upper bound on the black-hole fundamental (slowest damped) frequency 

min{Qti>} < ttTbh ■ (13) 

Thus the relaxation bound implies that a black hole must have (at least) one quasinormal resonance whose imaginary 
part conform to the upper bound (fiU)) . This mode would dominate the relaxation dynamics of the perturbed black 
hole and will determine its characteristic relaxation timescale. Taking cognizance of Eq. pip for the equatorial modes, 
and substituting n = for the fundamental resonance, one obtains min {Qui} = kTbh + 0(T 2 H ) [H M, 0, S, m. 



One therefore concludes that rapidly rotating Kerr black holes actually saturate the universal relaxation bound. 
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